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Based on renormalization group concepts and explicit mean field calculations we study the uni- 
versal contribution to the effective force and torque acting on an ellipsoidal colloidal particle which 
is dissolved in a critical fluid and is close to a homogeneous planar substrate. At the same closest 
distance between the substrate and the surface of the particle, the ellipsoidal particle prefers an 
orientation parallel to the substrate and the magnitude of the fluctuation induced force is larger 
than if the orientation of the particle is perpendicular to the substrate. The sign of the critical 
torque acting on the ellipsoidal particle depends on the type of boundary conditions for the order 
parameter at the particle and substrate surfaces, and on the pivot with respect to which the particle 
rotates. 



I. INTRODUCTION 

The confinement of the order parameter fluctuations 
in a critical fluid leads to an effective long-ranged inter- 
action between the confining walls and colloidal particles 
suspended in the fluid. The occurrence o f such a force 
was predicted by iFisher and de Gennesl in Il978l [l[ . It is 
called critical (or thermodynamic) Casimir force 0, 0], 
in analogy with the Casimir force [4| in quantum electro- 
dynamics where the force originates from the confined 
quantum fluctuations of the electromagnetic field Q. 

Since then the critical Casimir effect has attracted in- 
creasing theoretical and experimental attention. So far 
the theoretical investigations of the critical Casimir effect 
have been focused on the film geometry, realized either by 
homogeneous, planar and parallel walls (see, e.g., Refs. 
d 0, 1 ©, M, El M 0, EE MM and references 
therein) or by chemically patterned [181 ] or geometrically 
structured substrates pJI, as well as on spherical col- 
loidal particles (see, e.g., Refs. [2(1 HH, [H, [23| and refer- 
ences therein). Strong experimental evidences for critical 
Casimir forces have been obtained by studying wetting 
films near critical end points of quantum [2J, [2fl [26| or 
classical fluids [13, Ei| . 

It has been pointed out by IFisher and de Gennesl 
P, [29[ that such fluctuation-induced forces should lead 
to flocculation of colloidal particles if their solvent is a 
binary liquid mixture close to its consolute point. Such a 
solvent-mediated flocculation, which can be interpreted 
as indirect evidence for the critical Casimir force, has in- 
deed been observed for silica spheres suspended in a bi- 
nary liquid mixture of water and lutidine [30l . l3ll . l32l . [33| 
(see also Ref. [34| and references therein) as well as in 
other binary mixtures [H, [H,[37, 38]. However, only re- 
cently, the first direct experimental evidence of the criti- 
cal Casimir force has been reported concerning the force 
between a spherical colloidal particle and a homogeneous 
{39l ] or chemically patterned [lO, SI| wall. 

In the present study we analyze the critical Casimir 
effect for non-spherical colloidal particles. In this case, 
there is not only a force acting between particles or be- 



tween a particle and a wall but there is also a torque ex- 
erted on the particle. This may lead to interesting effects 
such as orientational ordering of non-spherical colloids in 
a critical solvent or anchoring of non-spherical particles 
at a wall. By varying the temperature the strength of 
this orientational interaction can be tuned and by chang- 
ing chemically the preferences of the surfaces for the two 
species forming the solvent one can choose the sign of the 
interaction [lg, 53, 511 . Motivated by this prospect we 
therefore extend our previous study of critical adsorption 
at a single non-spherical colloidal particle [12] to the case 
that in addition there is a planar wall present. 

This kind of orientation dependence of fluctuation- 
induced forces has recently been studied for quantum me- 
chanical Casimir forces [43|, EH . As far as torque due to 
critical fluctuations is concerned, in Ref. (4f| the critical 
Casimir torque on the confining walls of a wedge has been 
analyzed. Based on field-theoretic techniques the interac- 
tion of non-spherical particles, embedded into a solution 
of long polymers, with a planar wall has been investi- 
gated in the limiting case that the size of the particle is 
much smaller than the distance from the wall ("protein 
limit" ) and which, in turn, is assumed to be much smaller 
than the correlation length [H, 53, 51, 51] . These latter 
analyses show, in particular, that for a solution of ideal 
polymer chains the preferred orientation of the elongated 
colloidal particle changes from being perpendicular to be- 
ing parallel to the substrate surface upon decreasing the 
particle-wall distance [4(1 53], whereas in a solvent of 
self-avoiding chains the preferred orientation is the par- 
allel one for all distances [49| . 

The theoretical interest in the behavior of non- 
spherical colloidal particles is matched by an increased 
experimental interest, even with application perspec- 
tives (50j. Rodlike [Fl[ or disklike |52| architectures, 
dumbbell-shapedparticles [53|, HH , and particles with el- 
lipsoidal shape [53, [5(| have been synthesized and char- 
acterized. The size of these particles ranges between 10 
nm and 10 /im. Very recently the influence of an effective 
torque exerted by a non-critical one-component solvent 
on dumbbell-shaped particles has been revealed using de- 
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FIG. 1: Schematic side view of an ellipsoid near a planar 
wall. The ellipsoid semi-axes are R = Ri < R2 < ■ ■ ■ < Rd- 
Only the projection onto the (xi, xd) = {y,z) plane is shown. 
The angle between the long axis of the ellipsoid and the wall 
surface is denoted as a, and the closest distance between the 
wall and the surface of the ellipsoid is denoted as L, while 
the distance between the center of the ellipsoid and the wall 
is denoted as H . The full circle and the square denote two 
pivots which we use to calculate the torque (see Subsec. UlI C|l . 



polarized light scattering [57|. For this system, it has 
been demonstrated that the addition of small amounts of 
electrolyte has a significant impact on the rotational mo- 
tion and the aggregation stability of these non-spherical 
particles, while changes of the temperature typically re- 
sult in only minor changes of the effective interaction be- 
tween colloidal particles because this solvent is not close 
to a continuous phase transition. However, as stated 
above, near a critical point the effective interaction is 
expected to exhibit a very sensitive temperature depen- 
dence. 

The reminder of the paper is organized as follows. In 
Sec. |TT] we define the system under consideration and we 
introduce the scaling functions for the critical Casimir in- 
teractions. In order to calculate the force and the torque 
we use the stress tensor, as described in Subsec. IIII Al 
The results of the full mean field calculations are dis- 
cussed in Subsecs. llII Bl and IIII CI In Sec.|TV]we compare 
qualitatively the quantum-electrodynamic Casimir, crit- 
ical Casimir, and polymer depletion interactions acting 
on an ellipsoidal colloid close to a wall. Finally, in Sec.lVl 
we briefly summarize our results. 



II. MODEL 

In order to be able to discuss and to compare the be- 
havior of both non-spherical as well as spherical colloidal 
particles, we describe them in a unified way by consider- 



ing a particle with the shape of a hypercylinder: 
K. d ({Ri}) = jr = {xi,x 2 , ■ ■ ■ ,x D ) 
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where R — R\ < R 2 < ■ ■ ■ < Rd are the semi-axes of the 
hypercylinder (see Fig. [T]) . If R — Ri = ■ ■ ■ = R4 and 
d = D, the hypercylinder reduces to a hypersphere. If 
not all semi-axes are equal but d — D the hypercylinder 
is an ellipsoid. It is called a spheroid if only two semi- 
axes are different: in D — 3 one has a prolate spheroid 
(i?i = i?2 < -R3) and an oblate spheroid (i?i < R 2 = R3). 
In the case 1 < d < D, we have a spheroidal cylinder (or 
"spheroido-cylinder" ) , with a D — d dimensional hyper- 
axis. For reasons of simplicity, and in order to distinguish 
the cases of equal and different semi-axes, we shall use 
the notion of a disk (d = 2) and of a d-sphere (d > 3) 
if all finite axes are equal (note that they are cylinders 
for D > d), while the case in which some of the axes 
are different will simply be denoted as an ellipse or as 
a spheroid (or ellipsoid) depending on d (note that they 
are spheroido-cylindcrs in D > d). 

We consider the generalization of D to values differ- 
ent from three because the upper critical dimension (for 
the Ising universality class of the fluids discussed here) 
is D* = 4. We recall that the universal quantities such 
as critical exponents and universal scaling functions cal- 
culated within mean field theory (MFT) are exact in 
D > D* . (For D = D* one expects logarithmic correc- 
tions which we do not consider here.) Accordingly, our 
MFT results (see Sec. 1111 j) can be interpreted either as 
exact results in D > D* or as mean field approximations 
for the dimensions D = 3 or D = 2. 

The effective force / acting on a particle in the direc- 
tion normal to the substrate is the negative derivative of 
the free energy T of the confined fluid with respect to the 
closest distance L of the particle from the wall at fixed 
orientation. Close to a critical point it decomposes into 
the sum of an analytical background part and of a non- 
analytical (singular) part @, Q, [II (f° r further details see 
Ref. [HI]). According to the scaling behavior predicted 
by renormalization group theory the latter (divided by 
the "length" of the (D — d)-dimensional hyperaxis of the 
hypercylinder along which K,d is translational invariant), 
in the case of an ellipsoidal particle close to the wall, can 
be cast into the form (s = 1, • • • , d — 1) (22J: 



f(L,R 1 ,--- ,R d ,{a p },T) = 



J^D-d+l 

K ± (e ± = L/t ± ,A = L/R,{a p }; 

{S s = R S+1 /R}\, 



(2) 



where {a p } are d(d — l)/2 angles determining the ori- 
entation of the hypercylinder, £± = £g = |i| _!/ is the bulk 
correlation length in the disordered (+) and ordered (— ) 
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phase, respectively, t = (T — T c ) /T c with v as a standard 
bulk critical exponent, and £^ are non- universal ampli- 
tudes. L is the closest distance between the surface of the 
ellipsoid and the wall (see Fig. [1]), R = Ri, and K± are 
dimensionless universal scaling functions. / > (/ < 0) 
corresponds to repulsive (attractive) forces. 

In addition, a non-spherical particle experiences an ef- 
fective torque, the components of which are, in general, 
linear combinations of the (negative) derivatives of the 
free energy with respect to the angles {a p }. In the case 
shown in Fig. [TJ the torque has only one non-zero com- 
ponent describing the rotations within the yz plane; this 
component equals the negative derivative of the free en- 
ergy T with respect to the angle a. In general, the non- 
analytic contribution to the torque per "length" of the 
(D — e?)-dimensional hyperaxis of K,d takes on the scaling 
form 

t(L, Ri, ■ ■ ■ , R d , {a p }, T) = jpZi 

xM ± (e ± ,A,{a p };{5 s }), (3) 

where M± are the universal scaling functions of the 
torque. We note that t (and therefore M±) is an antisym- 
metric tensor of the second rank (the so-called 2-form) . 
In D dimensions t has, in general, (^2) = D\/(2(D— 2)!) 
components. However, because of the translational in- 
variance m (D - d) directions, only Q = d(d - l)/2 of 
them are non-zero (some of them may as well be zero 
depending on additional symmetries of ICa)- In D = 3 
dimension the conventional torque pseudo- vector t is ob- 
tained by taking the Hodge dual of t, i.e., t = *t; in 
Euclidian space t& = (1/2!) J2i j Uj^ijkt where e^k is the 
Levi-Civita symbol. 

In the literature two common choices are used to de- 
fine the distance between the particle and the wall: The 
surface-to-surface distance L (used, e.g., in Refs. [22H23I 
l4lj ]) and the distance H from the center of the particle 
to the wall (used, e.g., in Refs. [44]) |4y, [47|, [48|, |49|)> see 
Fig. HJ Correspondingly, one can also define two differ- 
ent pivots (points of rotation) with respect to which the 
particle rotates: The point on the surface of the particle 
closest to the wall (denoted by a circle in Fig. [1]) in the 
former case and the center of a particle (denoted by a 
square in Fig. [1]) in the latter case. While it is obvious 
that the force does not depend on the definition of the 
distance, the torque does depend on the choice of the 
pivot. It is easy to see, however, that the values of the 
torque corresponding to two different pivots, and i^ 2 \ 
are related by the simple equation (see, c.f., also Eq. (HH))) 

t (D = t ( 2 ) + r ( 12 ) A /, (4) 

where r^ 12 ) = — is the vector connecting the two 
pivot points located at and r^ 2 \ respectively. / is 
the vector of the force, and A denotes the wedge product 
(a multi-dimensional analogue of the three-dimensional 
cross product; its ij component is given by (r^ 12 ) A/)y = 



x\ fj — Xj fi, where x\ is the i-th component of 
r^ 12 ) . In the following we shall calculate the scaling func- 
tions K± and M± within MFT, and we will also compare 
and discuss in more detail the scaling functions for these 
two choices of the distance and of the pivot (see Sec. lIII Bl 
and[lHC|. 

III. MEAN-FIELD APPROXIMATION 

The standard Landau-Ginzburg-Wilson Hamiltonian 
for critical phenomena confined to a volume V is given 
by 

um = J v dv Qw) 2 + + ^ (5) 

augmented by boundary conditions [H, [5!| . In the case 
of a binary liquid mixture near its consolute point the 
order parameter <fi is proportional to the difference of the 
concentrations of its two species; r is proportional to t 
and u is a coupling constant which stabilizes H[(f>] in the 
ordered phase. For the critical adsorption fixed point 
[5!^ |. valid for confined fluids, the boundary conditions 
are 4> = +00 (or (fi — —00) at the surface of the col- 
loidal particle and at the wall, to which we refer to in 
the following as the "+" (or "— ") boundary condition. 
The semi-axes {Ri} of the hypercylinder (see Fig.[l| in- 
troduce additional length scales which might come into 
play via coupling constants of additional surface terms 
in the effective Hamiltonian [5!|. However, on the ba- 
sis of power counting one concludes that such terms are 
irrelevant at the ordinary transition (corresponding to 
the Dirichlet boundary conditions), where the surface en- 
hancement coupling of the term proportional to (fi 2 at the 
surface asymptotically dominates 1/Ri contributions to 
couplings of symmetry-preserving boundary terms. On 
the same footing we expect that the distinctive feature 
of the normal transition, i.e., the occurrence of symme- 
try breaking with a resulting asymptotic divergence of 
the order parameter at the surface is asymptotically not 
affected by curvatures; curvatures might perhaps influ- 
ence the cross-over between the ordinary and the normal 
transitions. Accordingly, we do not expect the aforemen- 
tioned asymptotic boundary conditions to be modified 
by curvature. 

Within MFT, the fluctuations of the order parameter 
<fi are neglected and only the order parameter configura- 
tion with the largest statistical weight, m = y/u/6{<fi), is 
taken into account. Minimization of Eq. |5| leads to the 
Euler-Lagrange equation 

Am = rm + to 3 . (6) 

Equation ([B]) is solved numerically as function of r us- 
ing the finite element method (see, e.g., Ref. [60]). We 
consider only hypercylinders K,d{{Ri\) with d = 2 and 
d = 3, for which the problem is effectively two- and three- 
dimensional, respectively (for /Q>4 one has a four and 
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higher dimensional problem which is difficult to solve nu- 
merically). For reasons of simplicity we also restrict our 
considerations to the case of only two different semi-axes, 
R = Ri = ■ ■ ■ = Rd-i < Rd, i-e., we consider a prolate 
spheroid in spatial dimension d (which is a spheroido- 
cy finder in spatial dimension D > d). We expect that the 
results for an oblate spheroid do not differ qualitatively. 
In Eq. ^ t can be expressed in terms of the bulk corre- 
lation length £± , which governs the exponential decay of 
the two-point correlation function in the bulk: r = 
for t > and r = -fl 2 / 2 for r < 0. 

We note that the contribution from the square gradi- 
ent term in Eq. © does not vanish in the case of the 
ellipsoid-wall geometry with symmetry-breaking bound- 
ary conditions. Therefore, one does not expect an addi- 
tional thermodynamic length [6f| to emerge and to affect 
finite-size scaling in dimensions D > D*. Such a thermo- 
dynamic length naturally emerges in the finite-size scal- 
ing analysis for D > D* if the standard mean-field theory 
in a finite volume exhibits an isolated zero mode which 
becomes "massless" at the bulk critical point due to the 
absence of a nonvanishing contribution from the square 
gradient term in Eq. ([5]). 



A. The stress tensor 

By using the stress tensor, the force and the torque 
can be calculated directly from the order parameter pro- 
file m(r). This has the advantage that one avoids the 
numerical difficulties of calculating differences of free en- 
ergies which attain large values due to the divergence of 
the order parameter profile at the surfaces. 

We consider an infinitesimal, local coordinate transfor- 
mation (k = 1, ■ ■ ■ , D): 



x'k = %k + X, Suj(r), 



(7) 



where the meaning of the index j depends on the kind 
of transformation and is specified below. The linear re- 
sponse of a system to such a coordinate transformation 
is 



J,k 



(8) 



where (d n m = dm/dx n ) 
d(d k m) 



Fig. [Q), divided by the "length" I = / d D ~ d x of the 
(D — G?)-dimensional hyperaxis of Kd is given by [23| 
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I dSH 

£ da 



1 



D 



- <h dS^2 %k ™k 

S k=l 



where within MFT T = min^ (%[</)]). S denotes the sur- 
face of Vo , h k is the fc-th component of its unit outward 
normal, dS = d D ~ 1 x, and Tjk is the conventional stress 
tensor: 



Tjk - 



S(d k m 



(11) 



One often adds to T the so-called "improvement" term 
(see, e.g., Refs. [H, Q 0, H), 



1-kl — 



ID -2 



4D- 



- [d k di-6 k iA]m 2 , 



(12) 



which ensures the scale and conformal invariance of the 
stress tensor at the critical point and renders it renormal- 
izable [62|, [H, [U [6^1 • However, the contributions from 
this term to 'observable quantities' like the force or the 
torque vanish. For instance, in the case of the force one 
obtains with the help of the Gauss-Ostrogradsky theo- 
rem (using here the summation convention) § s IiktiidS = 
J Vo diXikdV = because d{Ii k = 0. 

In order to calculate the torque, the coordinate trans- 
formation ([7]) is chosen to take the form of an infinitesi- 
mal rotation, i.e., 



Xk 



D 

E 

n=l 
n^k 



D 



X n Su! nk — X k + ^2 X k,nl SljJnU (13) 



i,n=l 
n<l 



where X k _ n i = xiSkn — x n 5 k i and the antisymmetry of 
5uj n k has been used. Note that here the index j in Eq. ([7]) 
consists of two indices, j — > (In), which denote the ro- 
tation plane. Now we choose 8u) yz (r) — a for r G Vo 
and 5u! n k(r) = otherwise (see Fig.[fJ. This renders the 
(only non-zero) component of the torque (per "length" 
of the (D — c?)-dimensional hyperaxis of /Q/): 



ty g i dsn _ 1 



k B T £ da £ 



j> dS ^2M kt y Z n k 

k=l 



H~da' 
(14) 



^2(d n m)Xj n CX jk . (9) w here the "angular momentum" tensor Mk, n l is 



In Eq. ((9]) C(m, dm) is the integrand of the Hamilto- 
nian H given by Eq. ([5]). In order to calculate the force, 
we use the coordinate transformation with Xj k = Sj k 
(j,k = 1, • • • ,D) and 5us z (r) = a for r £ Vq and 
6u>k(v) = otherwise, where Vq is a generalized hyper- 
cylinder enclosing /Q- Accordingly, the z-component of 
the force, which is the only non-zero component (see 
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k,nl 



Xl - X) 



T fc , 
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with Xf being the / component of the position vector of 
the pivot with respect to which the ellipsoid rotates. If we 



(15) 



choose x 



(i) 



,.< 2 ) 



+ xf 1 , where x^' is the I component 

(12) 

of the position vector of another pivot and x) is the I 



„(2) ! 



component of the vector r^ 12 ) = —r^ connecting the 
two pivots, then by using Eqs. (fTU|) and (jT5J) we obtain 



(1) _ ,(2) 



a 



,(12) 



fl 



,(12) 



fn 



(16) 



where /; is the I component of the force / (we recall that 
in our case only f z = / ^ 0), and denotes the (nl) 
component of the torque corresponding to the i-th pivot. 
Equation (fTo|) is the component version of Eq. (g]) . 



B. Casimir force 

We first study the influence of the anisotropy of the 
particle on the Casimir force for D — 4. In Fig. [2Kb) we 
compare the universal scaling functions K± , as functions 
of the rescaled distance 0± = L/£±, for a spheroidal 
particle JCd=d, and for two (d — 3)-spheres (see Eq. (p}) 
in the case of H — h boundary conditions which lead to 
attractive forces. In particular, we compare the two ex- 
treme orientations a = and a = tt/2 of a spheroid (i.e., 
the particle is oriented parallel and perpendicular to the 
wall, respectively, see Fig.[T]) with two (d = 3)-spheres: a 
small one with its radius equal to the smallest semi-axis 
(R = Ri), and a bigger one with its radius equal to the 
major semi-axis of the spheroid (R = Rd = 2Ri). It is 
interesting that the force is stronger for the small sphere 
((2) in Fig.[2{a)) than for the spheroid oriented perpen- 
dicular to the wall ((1) in Fig. Efa)). Indeed, within the 
Derjaguin approximation [66l [63] the force is inversely 
proportional to the square root of the Gaussian curva- 
ture G = k^k 1 - 2 ^ of the particle surface at the point 
closest to the wall (c.f., Eq. (fT7)) h are the principle 
curvatures. The square root of the Gaussian curvature 

1 /2 

of the sphere is G S p here = l/i?i while the principle cur- 
vatures of the spheroid at its elongated edge (a = tt/2) 



are k 



(i) 

spheroid 



spheroid 



R d /R\ so that G^Loidt" 

hero if R d > Rl 

(which is the present case, see Fig. [2{a)). For a = 
((3) in Fig. Eta)) one has G^LoidC" = V 2 ) = W 
which equals G^ ere = l/R d for (4) in Fig. Ufa). Thus 
the configurations (3) and (4) in Fig. H^a) have the same 
Gaussian curvature at the point closest to the wall so 
that within the Derjaguin approximation given by, c.f., 
Eq. (fT7|) both cases lead to the same critical Casimir 
force. Therefore the difference between the curves (3) 
and (4) in Fig. [2](b) highlights the shortcomings of the 
Derjaguin approximation (Eq. (|17[) i for this geometrical 
set-up. The scaling functions for intermediate orienta- 
tions and for H — h as well as H — boundary conditions are 
shown in Fig. [3] The minima of the scaling functions for 
H — h boundary conditions occur at T > T c . They move 
closer to T c upon changing the orientation of the particle 
from parallel to perpendicular. 

It is instructive to compare our full MFT results 
with the Derjaguin (or proximity force) approximation 
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FIG. 2: (color online) (a) Two spheres ((2) and (4)) and one 
spheroid with perpendicular (a = tt/2; (1)) and parallel (a = 
0; (3)) orientation at the same surface-to-surface distance L 
from the wall, (b) The universal scaling functions K± — 
K±{e±,A = 2, a; S = 2)/\K+{@+ = 0, A = 2; S = 1)| as a 
function of the rescaled reduced temperature 0± = L/£± = 
\t\"L/(£, where 8 = Rd/R- The curves (2) and (4) correspond 
to (d — 3)-spheres with radii R — Ri — L/2 and R — Rd = 
2Ri = L, respectively. The curves (1) and (3) correspond to 
a prolate spheroid with R = R\ = R% < Rd=3 (spheroido- 
cylinder /C3 in spatial dimension D > 4) for A = L/R = 2, 
S — Rd=s/R = 2, and the orientations a = tt/2 and a = 0, 
respectively, where a is the angle between the main axis of the 
spheroid and the wall (see Fig. [1} . In all cases ++ boundary 
conditions are imposed. The force is expressed in units of the 
absolute value of the Casimir force for the sphere with radius 
R = L/2 at the critical point (0+ = 0) and for ++ boundary 
conditions (accordingly, at 0+ = line (2) is normalized to 
-1). 



[66l |67| . which is applicable for large particles close to 
the wall. In lowest order in curvatures one obtains 

\(d-l)/2 _ 

K±(e,A,a; S) = G~ 1/2 (A, 8, a) 



r((d + l)/2) 



( u - l)(<*- 3 )/ a du, (17) 



where K± is the scaling function for the plate-plate ge- 
ometry, r(a;) is the Gamma function, and G(A, S, a) = 
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(a) 



FIG. 3: (color online) The universal scaling function K+ — 
K+(Q+,A = 2, a; S = 2)/\K+{&+ = 0, A = 2, a = 0; 8 = 
2)| as a function of the rescaled reduced temperature 0+ = 
L/£+ = i" L/£,o f° r a prolate spheroid R = R\ = R2 < Rs 
(spheroido-cylinder /C<j=3 in D > d) for A — L/R — 2 and 
5 = R3/R = 2, and for three values of the angle a between 
the main axis of the spheroid and the wall (see Fig. [1} . The 
positive (negative) curves correspond to +— (H — h) boundary 
conditions on the spheroid and the wall leading to repulsion 
(attraction). The force is expressed in units of the absolute 
value of the Casimir force for a = at the critical point (0+ = 
0) and for ++ boundary conditions. Accordingly, at 0+ = 
the bottom curve is normalized to —1. (Note that here the 
normalization differs from the one used in Fig. [5]) The dots 
denote the minima of the scaling function. Dotted and dash- 
double dotted lines show the Derjaguin approximation (DA) 
for a — and a = tt/2, respectively, and for +— boundary 
conditions. According to Eqs. (|17|l and (|18ft these two curves 
differ only by an overall scale factor S 2 — 4. 



j_d-\ FJ^ = * (where are the principle curvatures) 
is the dimensionless Gausian curvature of the particle 
surface at the point closest to the wall. In the case of the 
spheroid /Q=3 one has 



G(A,6,<p(a)) 



A 2 S 



2x2 



(1 + (5 2 - 1) cos 2 <p(a)Y 



(18) 



where tp is the parametric latitude denoting the 
position on the spheroid surface parameterized as 
(i?i cos ip cos tp, R\ cos sin '0, Rdsimp) with ip as the 
longitude. The value of ip at the point closest to the 
wall depends on a; ip{a = 0) = and tp(a = tt/2) = tt/2. 

The results within the Derjaguin approximation (us- 
ing the mean field scaling function K±^) are shown in 
Fig. [3] by dotted and dash-double dotted lines. In or- 
der to avoid overloading Fig. [3J the Derjaguin results are 
presented only for two angles, a = and a — tt/2. For 
these two angles the Derjaguin results differ only by an 
overall scale factor S 2 with S 2 = 4 here. For both angles 
the discrepancy with the full mean field results is rather 
large. The reason is that in the case studied here the 
ratio A = L/R = 2, while the Derjaguin approximation 
is supposed to be valid in the limit of small distances 
(or big particles), i.e., for A — ► 0. It is interesting to 




- a = 0, L/R = 2 

- a = tt/2, L/R = 2 
a = 0, L/R = 3 




FIG. 4: (color online) (a) Comparison between three ellipse 
configurations such that the configurations (1) and (2) have 
the same closest distance L = 2R, while configurations (2) 
and (3) have the same distance H = 4R between the center 
of the particle and the wall. For the configuration (1) H = 3R 
and for the configuration (3) L — 3R. (b) The universal scal- 
ing function K+ = K+(G+,A,a; S = 2)/\K+{0+ = 0, A = 
2, a — 0; 6 — 2)| as a function of the reduced temperature 
Q+ — L/£+ for an ellipse (d = 2). The dashed lines (1) rep- 
resent the results for an ellipse oriented parallel to the wall 
(a = 0, L/R = 2, H/R — 3) at the same closest distance 
between its surface and the wall as an ellipse which is ori- 
ented perpendicular to the wall (a — tt/2, L/R = 2, H/R — 4, 
solid lines (2)). The dash-dotted lines (3) represent the re- 
sults for an ellipse oriented parallel to the wall (a — 0, 
L/R — 3, H/R = 4) and at the same distance between its 
center of mass and the wall as compared to the ellipse which 
is oriented perpendicular to the wall (solid lines (2)). The 
positive (negative) curves correspond to +— (++) boundary 
conditions on the ellipse and at the wall. The force is ex- 
pressed in units of the absolute value of the Casimir force for 
a — and A = L/i? = 2at the critical point (O+ = 0) and 
for ++ boundary conditions. 



note that the Derjaguin approximation underestimates 
the strength of the force for a = tt/2 and overestimates 
it for a — 0. 

For the comparison of the Casimir force acting on el- 
lipsoids with various orientations and at a fixed distance 
from the wall in ^-direction, there are two interesting 
choices for this distance: The closest distance L between 
the wall and the surface of the ellipsoid, and the distance 
H between the center of the ellipsoid and the wall (see 
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Fig. [T]) . Figure 0|b) displays the corresponding scaling 
function for an ellipse (d = 2) which is oriented perpen- 
dicular to the wall (a = n/2,L/R = 2, H/R = 4, solid 
lines; (2) in Fig. 01 together with the scaling function 
for an ellipse which is oriented parallel to the wall at the 
same distance L (a = 0, L/R = 2, H/R = 3, dashed 
lines; (1) in Fig. 2]) and the scaling function for an ellipse 
which is oriented parallel to the wall at the same distance 
H (a = 0, L/R = 3, H/R = 4, dash-dotted lines; (3) in 
Fig. rj|. If the ellipse is orientated perpendicular to the 
wall (solid lines), the magnitude of the Casimir force is 
smaller than if the ellipse is parallel to the wall at the 
same closest distance between the wall and the surface of 
the ellipse (dashed lines). Compared to an ellipse which is 
parallel to the wall at the same distance between the cen- 
ter of the ellipse and the wall (dash-dotted lines) the mag- 
nitude of the Casimir force acting on the perpendicular 
oriented ellipse is larger (solid lines). We have confirmed 
numerically that these relations between the strengths of 
the fluctuation induced forces are valid for arbitrary val- 
ues of A — L/R including the limits A — > and A>1. 
In the latter limit these results are in agreement with the 
predictions of a small particle operator expansion (68| 
which can be used near the critical point (0 + — > 0) as 
long as all semi-axes of the ellipse are much smaller than 
all other lengths such as the correlation length or the 
distance between the surface of the particle and the wall. 

In Fig. [5] the scaling function K + with a suitable nor- 
malization is plotted as a function of the angle a between 
the main axis of the spheroid and the wall (see Fig. [T]) 
for a fixed rescaled minimal distance = const. 

As expected, the force is larger for more elongated col- 
loids (Rd/Ri large) but only for small a (i.e., if the 
colloids are almost parallel to the wall). If the col- 
loids are tilted more towards the perpendicular orien- 
tation (a > 30°), the opposite trend is observed. As 
discussed at the beginning of this subsection, this is due 
to the shapes of the prolate spheroids: for the perpen- 
dicular orientation G~V5 oid (a = 7r/2) = R 2 / Rd for 



* and Gj/ p l c (7r/2) 



R/R l J 2 for d = 2, which 
decrease upon increasing Rd, whereas for the parallel 



orientation one has G 

-1/2 



-1/2 
spheroid 



(o) 



Rd for d = 3 and 



^ellipse (0) = Rd/R 1 ^ 2 for d = 2 which increase upon in- 
creasing Rd- The scaling functions for the hypercylinders 
^3 (Fig-IS[a)) and K,2 (Fig. [Sfb)) are rather similar. 



C. Casimir torque 

The torque scaling function M + (see Eqs. ([3]), (fl4|) . and 
(JT5J)) as a function of the angle a (Fig. [I} between the 
main axis of an elongated colloid and the wall is presented 
in Fig. [S] for both ++ and H — boundary conditions and 
for a fixed rescaled minimal distance L/£+ = const with 
the pivot point (see Eq. (|15p) taken to be the point 
closest to the wall (denoted as a circle in Fig. [T]). The 
scaling function is positive for H — h and negative for H — 



1 1 1 1 1 1 1— 

\. a = 3, A = 2, <d + = 2 


(a) 


\ R 'Ft - 


- O 


■ n. _ \ Rd/ R = 


= 3 


v \ Rd/R - 


= 2 






1 1 





tt/6 tt/3 tt/2 

a 



d = 2, A = 1, 9+ = 3 




Rd/R = 5 

Rd/ R = 3 

Rd/R = 2 



tt/6 



Q 



tt/3 



tt/2 



FIG. 5: (color online) The universal scaling function K+ = 
K+(e+,A,a; 8)/K+(Q+ = 0, A; 5 = 1) for +- boundary 
conditions as a function of the angle a between the main axis 
of the spheroid and the wall (see Fig. fl) for three values of 
the ratio 5 = Rd/R- In (a) d = 3, 0+ = = 2, and 

A = L/R = 2, while in (b) d = 2, 6+ = 3, and A = 1. The 
force is expressed in terms of the Casimir force at the critical 
point for +— boundary conditions for a sphere (d — 3) with 
radius R = Ri — L/2 in (a) and for a disk (d — 2) with radius 
R = L in (b). All curves correspond to keeping temperature 
and the minimal distance L fixed upon varying the orientation 
(see Fig. [[}. 



boundary conditions. This means that in the ++ case the 
configuration is optimal (i.e., the free energy is lowest) 
if the colloid is parallel (a — 0) while in the H — case 
if it is perpendicular to the wall (a — tt/2). As one 
may expect, the torque vanishes for a — and a = 
7r/2, which for H — (H — (-) boundary conditions correspond 
to a maximum (minimum) and a minimum (maximum) 
of the free energy, respectively. Also as expected, the 
magnitude of the torque increases upon increasing the 
ratio Rd/Ri (for fixed R\/L). 

As we have already mentioned, one can introduce dif- 
ferent pivots with respect to which the torque is exerted 
on the ellipsoid, such as the point of closest approach to 
the wall (denoted as a circle in Fig. [T]) and the center of 
the ellipsoid (denoted as a square in Fig. The latter 
is more convenient to use if the ellipsoid is far from the 
substrate. In this case the motion of the ellipsoid can be 
described in terms of its center of mass (which we assume 
to coincide with the geometrical center); accordingly the 
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FIG. 6: (color online) The universal torque scaling function 
M+ = M + (Q + ,A,a; S)/\K+(Q+ = 0, A; 5 = l)j as a func- 
tion of the angle a between the wall and the main axis of 
an elongated colloid (Fig. [1} for three values of the ratio 
S = Rd/R. In (a) d = 3, A = L/R = 2, and 6+ = L/£ + = 2, 
while in (b) d = 2, A = 1, and <9 + = 3. The positive (neg- 
ative) curves correspond to +— (++) boundary conditions. 
M+ > (M+ < 0) implies that the torque acts as to increase 
(decrease) a. The torque is expressed in terms of the absolute 
value of the Casimir force for ++ boundary conditions at the 
critical point for a (d = 3)-sphere with radius R — Ri — L/2 
in (a) and for a disk (d — 2) with radius R — L in (b). The ex- 
trema are marked by full dots; their positions shift to smaller 
values of a upon increasing Rd/R- All curves correspond to 
keeping temperature and the minimal distance L fixed upon 
varying the orientation, and the pivot is taken to be the point 
closest to the wall (denoted as a circle in Fig. [JJ. 




2 4 6 8 10 



t ■ 1 ■ 1 ■ r 




2 4 6 8 10 

e+ = L/e+ 

FIG. 7: (color online) The universal torque scaling func- 
tion M+ = A4 + (Q+,A,a; S = 2)/K { } ) (@ + = 0) as a func- 
tion of temperature for a spheroid (spheroido-cylinder !Cd=3 
in D > d) with 5 — Rd/R = 2 for ++ boundary conditions 
and for two values of the angle a between the main axis of the 
spheroid and the wall (Fig. [TJ. In (a) the surface-to-surface 
distance L is fixed with A = L/R = 2 and the torque is calcu- 
lated with respect to the point of closest approach (denoted 
as a circle in Fig. [JJ . In (b) the distance H between the wall 
and the center of the ellipsoid is fixed with H/R = 4 and the 
torque is calculated with respect to the center of the ellipsoid 
(denoted as a square in Fig.[TJ. The insets show M+ as a func- 
tion of a for L/£+ = 2 = const in (a) and H/£+ = 2 = const 
in (b). The symbols in the insets correspond to the values 
of the angle a represented in the main plots. The torque is 
expressed in term of the absolute value of the Casimir force 
for ++ boundary conditions at the critical point for the film 
geometry. In both (a) and (b) the torque is strongest at a 
distinct temperature above (and not at) T c . 



orientational degrees of freedom of the ellipsoid should 
consistently be described also with respect to the center 
of the ellipsoid. If the ellipsoid is sufficiently close to the 
substrate, it is more convenient to monitor the orienta- 
tions of the ellipsoid with respect to the point of closest 
approach to the wall. One reason is that in the region 
close to the substrate not all orientations of the ellipsoid 
with respect to its center of mass are allowed because 
the particle cannot penetrate the substrate (c.f., Sec. IIVI 
and Fig. |UJa)). One can also imagine that the ellipsoid 
is trapped in a relatively shallow potential well of optical 
tweezers. In this situation the rotation of the ellipsoidal 
particle occurs naturally with respect to the point closest 
to the wall. 



In Figs. Ul a) and[7Jb) the torque scaling function M + 
with a suitable normalization is plotted for these two 
choices of the pivot for a few values of the angle a be- 
tween the main axis of the ellipsoid and the wall and 
for H — h boundary conditions. In Fig. [T^a) the torque is 
calculated with respect to the point closest to the wall, 
and consistently we consider a constant minimal rescaled 
distance L/£+ = const for all angles. In Fig. UH/b) the 
torque is calculated with respect to the center of the el- 
lipsoid, so that the rescaled distance between the center 
of the ellipsoid and the wall is fixed, i.e., H/^ + = const. 
Note that the scaling function is negative in the former 
and positive in the latter case. This means that if the 
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surface-to-surface distance L is kept fixed (for instance, 
by optical tweezers or by the wall), the optimal configu- 
ration of an ellipsoidal colloid is to be parallel to the wall 
(a = 0). However, if the ellipsoid rotates with respect 
to its center of mass, the optimal configuration is to be 
perpendicular to the wall (a — 7r/2). In the H — case the 
situation is reverse. 



IV. COMPARISON WITH THE 
QUANTUM-ELECTRODYNAMIC CASIMIR 
INTERACTION AND THE POLYMER INDUCED 
DEPLETION INTERACTION 

The quantum-electrodynamic Casimir interaction 43l 
9 and the polymer induced depletion interaction |46l . 
471 . |48L |49( lead to pronounced effects on the orientational 
ordering of non-spherical particles, too. The schematic 
presentation in Fig. [5] illustrates the influence of these 
two interactions and of the critical Casimir interaction 
on the orientation of a prolate ellipsoidal particle near a 
planar wall. For the purpose of the discussion, the region 
adjacent to the wall can be divided into three sub-regions 
denoted as I, II, and III (see Fig. [8J) . In region I not all 
orientations of the prolate ellipsoid are allowed because 
the particle cannot penetrate the wall. In the case of 
an attractive critical Casimir interaction the prolate el- 
lipsoid reaches its most favorable configuration of lying 
parallel to the wall by tilting such that for a prescribed 
distance of its center from the wall its optimum angle 
in region I is the one for which it is in touch with the 
wall (Fig. [5Ja)). The influence of both the quantum- 
electrodynamic Casimir interaction and the polymer in- 
duced depletion interaction on the orientation of an el- 
lipsoid in region I has not yet been studied. For larger 
distances from the wall, in regions II and III, the critical 
Casimir torque drives the prolate ellipsoid into an orien- 
tation perpendicular to the wall (Fig.[5Ja)). In the case of 
the quantum-electrodynamic Casimir interaction and of 
the polymer induced depletion interaction the preferred 
orientation in re gion II is perpendicular [44[ (Fig. [Hb)) 
and parallel [MEl EU (Fig. [gfc)) to the wall, respec- 
tively. Upon further increasing the distance from the wall 
a change of the preferred orientation has been found in 
the case of a scalar model of quantum-electrodynamics 
with Neumann boundary conditions both at the surface 
of the particle and at the wall [44(| (Fig. Hfb)) and in 
the case of ideal polymers acting as depletion agents [46[ 
(Fig. He)). 

The polymer induced depletion interaction has been 
calculated [4(| [48|, [49j in the so-called protein limit in 
which the size of the ellipsoid is small compared to the 
polymer size characterized by the radius of gyration. 
By using the small particle operator expansion, it has 
been shown that the extension of region II along the z- 
direction is given by the radius of gyration of the poly- 
mers acting as depletion agents [III [43, [H, H§| . However, 
the presently available non-spherical colloidal particles 



are larger than typical polymers and monitoring small 
particles by optical techniques is very difficult. There- 
fore, it would be rewarding to study theoretically and 
experimentally the polymer induced depletion interaction 
beyond the small particle limit. 

We emphasize that the preferred orientations of a pro- 
late ellipsoid due to the quantum-electrodynamic Casimir 
interaction (Fig. [5£b)) have been obtained by consider- 
ing a scalar model instead of the actual vectorial elec- 
tromagnetism. In the actual case of the electromagnetic 
field, which implies Dirichlet boundary conditions, and 
a slightly deformed sphere, the preferred orientation of 
a prolate ellipsoid is the one perpendicular to the wall 
in both regions II and III [44|. Future work in this 
area may focus on the understanding of the influence of 
the quantum-electrodynamic Casimir interaction on el- 
lipsoids of arbitrary eccentricities and at small distances 
from the wall. 



V. CONCLUSIONS 

We have investigated the critical Casimir effect for sin- 
gle non-spherical colloidal particles immersed in a fluid 
near its critical point and exposed to a laterally homoge- 
neous planar wall. For an ellipsoidal colloidal particle the 
resulting critical Casimir force and torque can be charac- 
terized by universal scaling functions K± (Eq. ((T|)) and 
M± (Eq. ©), respectively, which depend on the dimen- 
sionless scaling variables <d± — L/^±, A = L/Ri, the 
angle a between the main axis of the ellipsoid and the 
wall, and the ratios <5j_i = R4/R1, where R\ is the small- 
est semi-axis of the ellipsoid, Ri (i = 2, • • ■ , d) are d — 1 
remaining semi-axes, L is the closest distance between 
the surface of the ellipsoid and the wall (see Fig. [T]) , and 
£± is the bulk correlation length above (+) and below (— ) 
the critical point. The scaling functions have been cal- 
culated within mean-field theory, which represents the 
leading order term in a systematic e = 4 — D expan- 
sion. The dependence of the scaling functions K± on the 
scaling variable Q± exhibits behaviors which are qualita- 
tively similar to those of a spherical colloidal particle (see 
Figs. [2] and [3]). The strength of the force depends on the 
orientation of the colloidal particle relative to the wall: 
for elongated colloids it is stronger if the colloid is ori- 
ented parallel to the wall at the same surface-to-surface 
distance L (see Figs. 2] and [3]). If ellipsoidal colloids are 
oriented perpendicular to the wall, the force is stronger 
for shorter colloids (see Fig. [5]). We note, however, that 
the latter effect is due to the specific shape of the ellip- 
soid at its elongated edge (see Subsec. IIIIB|) . while the 
former effect is more general. 

The sign of the universal torque scaling function M + 
depends on the boundary conditions and on the pivot 
with respect to which the particle rotates. Thus, if the 
pivot is chosen to be the point on the particle surface 
closest to the wall (denoted as a circle in Fig. [T]), the 
scaling function is positive for equal (++) and negative 
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FIG. 8: (color online) Illustration of the preferred orientation of a prolate ellipsoid at a fixed distance of its center from a 
planar wall in the z-direction. Only the projection onto the plane of the figure is shown. The centers of the ellipsoids are 
chosen to lie within the plane of the figure, (a) In the case of the attractive critical Casimir interaction the ellipsoid is oriented 
perpendicular to the wall in regions II and III, while it prefers to touch the wall in region I because in that configuration the 
corresponding Casimir force is stronger than for other orientations, (b) In the case of a scalar model of quantum-electrodynamics 
and for Neumann boundary conditions on the surface of the ellipsoid, a change of the preferred orientation has been found 
for Neumann boundary conditions at the wall [44J. For Neumann boundary conditions on the particle surface and Dirichlet 
boundary conditions at the wall the energy is minimal if the ellipsoid is perpendicular to the wall in both regions II and III. 
In the case of mixed Dirichlet and Neumann boundary conditions the quantum-electrodynamic Casimir force is repulsive, (c) 
Using the small particle operator expansion it has been shown that the attractive polymer induced depletion interaction leads 
to a change of the preferred orientation of the ellipsoid in the case of the ideal pol ym ers, while the favorable orientation is 
the parallel one in both regions II and III in the case of self-avoiding polymers [46I . 14a . |49I |. The preferred orientation of an 
ellipsoidal particle in region I has not yet been studied for the quantum-electrodynamic Casimir force or the polymer depletion 
induced interactions. The extension of region II in the ^-direction is of the order of the length of the long axis of the ellipsoid 
in (b) and given by the radius of gyration of the polymers in (c). 



for opposing (-1 — ) boundary conditions at the wall and 
at the particle surface (see Figs. [6] and Fig. 0(a)). This 
means that at the same closest distance L between the 
surfaces of the colloid and the wall an elongated particle 
tends to orient itself parallel to the wall for equal bound- 
ary conditions and perpendicular to the wall for opposing 
boundary conditions. We expect this situation to be re- 
alized if a particle is close to the wall or, for instance, 
if it is trapped by optical tweezers. An opposite effect 
is observed if the center of the particle is kept fixed (see 
Fig. [7]Jb)). In this case an elongated colloid prefers an 
orientation perpendicular to the wall for equal boundary 
conditions and parallel to the wall for opposing boundary 



conditions. It is worthwhile to note that this conclusion 
agrees qualitatively with the asymptotic results obtained 
from the small-particle operator expansion [68[ (see also 
Pig- Ha)). 
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